PARTIAL REGULARITY FOR OPTIMAL TRANSPORT MAPS 



GUIDO DE PHILIPPIS AND ALESSIO FIGALLI 

Abstract. We prove that, for general cost functions on R n , or for the cost d 2 /2 on a Riemannian 
manifold, optimal transport maps between smooth densities are always smooth outside a closed 
singular set of measure zero. 



1. Introduction 

A natural and important issue in optimal transport theory is the regularity of optimal transport 
maps. Indeed, apart from being a typical PDE/analysis question, knowing whether optimal maps 
are smooth or not is an important step towards a qualitative understanding of them. 

It is by now well known that, for the smoothness of optimal maps, conditions on both the cost 
function and on the geometry of the supports of the measures are needed. 

In the special case c(x, y) = \x — y\ 2 /2 on W 1 , Caffarelli [3], SI E] proved regularity of optimal 
maps under suitable assumptions on the densities and on the geometry of their support. More 
precisely, in its simplest form, Caffareili's result states as follows: 

Theorem 1.1. Let f and g be smooth probability densities, respectively bounded away from zero 
and infinity on two bounded open sets X and Y , and let T : X — >■ Y denote the unique optimal 
transport map from f to g for the quadratic cost \x — y\ 2 /2. IfY is convex, then T is smooth inside 
X. On the other hand, ifY is not convex, then there exist smooth densities f and g (both bounded 
away from zero and infinity on X and Y , respectively) for which the map T is not continuous. 

A natural question which arises from the previous result is whether one may prove some partial 
regularity on T when the convexity assumption on Y is removed. In [16|. [18] the authors proved 
the following result: 

Theorem 1.2. Let f and g be smooth probability densities, respectively bounded away from zero 
and infinity on two bounded open sets X and Y , and let T : X — >■ Y denote the unique optimal 
transport map from f to g for the quadratic cost \x — y\ 2 /2. Then there exist two open sets X' C X 
and Y' C Y , with \X \ X'\ = \Y \ Y'\ = 0, such that T : X' —>Y' is a smooth diffeomorphism. 

In the case of general cost functions on M n , or when c(x,y) = d(x,y) 2 /2 on a Riemannian 
manifold M {d{x,y) being the Riemannian distance), the situation is much more complicated. 
Indeed, as shown by Ma, Trudinger, and Wang [33], and Loeper [31], in addition to suitable 
convexity assumptions on the support of the target density (or on the cut locus of the manifold 
when supp(g) = M [24]), a very strong structural condition on the cost function, the so-called 
MTW condition, is needed to ensure the smoothness of the map. 

More precisely, if the MTW condition holds (together with some suitable convexity assumptions 
on the target domain), then the optimal map is smooth |35 [ 1361 121] 130 } I19j. On the other hand, if 
the MTW condition fails at one point, then one can construct smooth densities (both supported 
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on domains which satisfy the needed convexity assumptions) for which the optimal transport map 
is not continuous [3l] (see also [15]). 

In the case of Riemannian manifolds, the MTW condition for c = d 2 /2 is very restrictive: indeed, 
as shown by Loeper [31], it implies that M has non-negative sectional curvature, and actually it is 
much stronger than the latter [28} 123]. In particular, if M has negative sectional curvature, then 
the MTW condition fails at every point. Let us also mention that, up to now, the MTW condition 
is known to be satisfied only for very special classes of Riemannian manifolds, such as spheres, their 
products, their quotients and submersions, and their perturbations [32], [22j [TOj, [29J EU [20l E] , and 
for instance it is known to fail on sufficiently flat ellipsoids [23] . 

The goal of the present paper is to show that, even without any condition on the cost function 
or on the supports of the densities, optimal transport maps are always smooth outside a closed 
singular set of measure zero. In order to state our results, we first have to introduce some basic 
assumptions on the cost functions which are needed to ensure existence and uniqueness of optimal 
maps. As before, X and Y denote two open subsets of W l . 

(CO) The cost function c : X x Y — > R is of class C 2 with ||c||c<2pf X y) < oo. 
(CI) For any x £ X, the map Y B y \— > —D x c(x,y) G R n is injective. 
(C2) For any y EY, the map —D y c(x, y) G R n is injective. 

(C3) det(D xy c)(x,y) / for all (x,y) G X x Y. 

Here are our main results: 



Theorem 1.3. Let X, Y C W 1 be two bounded open sets, and let f : X — > R + and g : Y — >• R + 
be two continuous probability densities, respectively bounded away from zero and infinity on X 
and Y. Assume that the cost c : X x Y — > R satisfies (C0)-(C3), and denote by T : X — > Y 
the unique optimal transport map sending f onto g. Then there exist two relatively closed sets 
Tix C X, Sy C Y of measure zero such that T : X \ — > Y \ Xy is a homeomorphism of class 
C°£ for any f3 < 1. In addition, if c G C^ 2 ' a (X xY), f G C^(X), and g G C^{Y) for some 
k > and a G (0, 1), then T : X \ T,x — > Y \ Sy is a diffeomorphism of class C^ 1,a . 

Theorem 1.4. Let M be a smooth Riemannian manifold, and let f,g : M — >■ R + be two continuous 
probability densities, locally bounded away from zero and infinity on M . Let T : M —> M denote 
the optimal transport map for the cost c = d 2 /2 sending f onto g. Then there exist two closed sets 
Sx, Sy C M of measure zero such that T : M \ T,x — >■ M \ Sy is a homeomorphism of class Cj^T 
for any f3 < 1. In addition, if both f and g are of class C k,a , then T : M \ T^x — > M \ Sy is a 
diffeomorphism of class C^ 1 '". 

The paper is structured as follows: in the next section we introduce some notation and prelim- 
inary results. Then, in Section [31 we show how both Theorem 11.31 and Theorem 11.41 are a direct 
consequence of some local regularity results around differentiability points of T, see Theorems 14.31 
and 15.31 Finally, Sections H] and [5] are devoted to the proof of these local results. 
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2. Notation and preliminary results 

Through a well established procedure, maps that solve optimal transport problems derive from 
a c-convex potential, itself solution to a Monge- Ampere type equation. 

More precisely, given a cost function c : X x Y — > R, a function u : X — > R is said c-convex if it 
can be written as 

(2.1) u(x) = sup {-c(x,y) + X y } , 

y eY 

for some constants A 9 G MU {— oo}. 

Similarly to the subdifferential for convex function, for c-convex functions one can talk about 
their c-subdifferential: if u : X — > R is a c-convex function as above, the c- subdifferential of u at x 
is the (nonempty) set 

(2.2) d c u{x) :={yeY: u(z) > -c(z,y) + c(x,y) + u(x) Vz£l}. 
If xo E X and yo £ 9 c it(a;o), we will say that the function 

(2.3) C XOi y (-) := -c(-,y ) + c(x ,y ) + u(x ) 

is a c-support for it at xo- We also define the Frechet subdifferential of u at x as 
<9~-u(x) := {p £ R" : u(z) > u(x) + p ■ (z - x) + o(\z - x\)} . 
We will use the following notation: if E C X then 

d c u(E) := (J c u(sc), a _ u(£:) := J ^n(x). 

It is easy to check that, if c is of class C , then the following inclusion holds: 

(2.4) y e d c u(x) => —D x c(x,y) E d~u(x). 

In addition, if c satisfies (C0)-(C2), then we can define the c-exponential map: 

(2.5) foranyxGJ^Gy.pGl", f c-ex Px .(p) = y O p = -Dc(x, y ) 

v ; J ' w ' ^ ' | c*-exp y (p) = x <^ p = -DyC{x,y) 



Using (|2.5p . we can rewrite (|2.4p as 

(2.6) d c u(x) C c-exp x (5~w(x)) . 

Notice that, if c E C 1 and y is bounded, it follows immediately from (|2.ip that c-convex functions 
are Lipschitz, so in particular they are differentiable a.e. 

The following notation will be convenient: given a c-convex function u : X — > R, we define (at 
almost every point) the map T u : X — > Y as 

(2.7) T u {x) :=c-exp x (Vu(x)). 

(Of course T u depends also on c, but to keep the notation lighter we prefer not to make this 
dependence explicit. The reader should keep in mind that, whenever we write T u , the cost c is 
always the one for which u is c-convex.) 

Finally, let us observe that if c satisfies (CO) and Y is bounded, then it follows from (|2.ip that u 
is semiconvex (i.e., there exists a constant C > such that u + C\x\ 2 /2 is convex, see for instance 
|13j). In particular, by Alexandrov's Theorem, c-convex functions are twice differentiable a.e. (see 
\37\ Theorem 14.25] for a list of different equivalent definitions of this notion). 

The following is a basic result in optimal transport theory (see for instance [371 Chapter 10]): 
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Theorem 2.1. Let c : X x Y — > R satisfy (CO) -(CI). Given two probability densities f and 
g supported on X and Y respectively, there exists a c- convex function u : X — > R such that 
T u : X — 7- y is the unique optimal transport map sending f onto g. 

In the particular case c(x,y) = —x • y (which is equivalent to the quadratic cost \x — y\ 2 /2), 
c-convex functions are convex and the above result takes the following simple form [2]: 

Theorem 2.2. Let c(x,y) = —x ■ y. Given two probability densities f and g supported on X and 
Y respectively, there exists a convex function v : X — > R such that T v = V'u : X — > Y is the unique 
optimal transport map sending f onto g. 

Although on Riemannian manifolds the cost function c = d 2 /2 is not smooth everywhere, one can 
still prove existence of optimal maps |34L I13|, [T7] (let us remark that, in this case, the c-exponential 
map coincides with the classical exponential map in Riemannian geometry): 

Theorem 2.3. Let M be a smooth Riemannian manifold, and c = d 2 /2. Given two probability 
densities f and g supported on M, there exists a c-convex function u:M->IU {+00} such that 
u is differentiable f -a.e., and T u (x) = exp x (Vu(x)) is the unique optimal transport map sending f 
onto g. 

We conclude this section by recalling that c-convex functions arising in optimal transport prob- 
lems solve a Monge- Ampere type equation almost everywhere, referring to [H Section 6.2], [3T|, 
Chapters 11 and 12], and [15] for more details. 

Whenever c satisfies (C0)-(C3), then the transport condition (T u )j/ = g gives 

(2.8) |det(ZXTjx))| = -M- a.e. 

g(T u (x)) 

In addition, the c-convexity of u implies that, at every point x where u is twice differentiable, 

(2.9) D 2 u(x) + D xx c(x, c-exp^, (Vit(x))) > 0. 
Hence, writing (12.7P as 

-D x c(x,T u (x)) = Vu(x), 
differentiating the above relation with respect to x, and using (|2.8|) and (|2.9|) . we obtain 
(2.10) 



det(D 2 u(x) + D xx c(x,c-exp x (Vu(x))) J = | det (D xy c(x,c-ex.p x (\7u(x)))) 



g(c-exp x (Vu(x))) 

at every point x where u it is twice differentiable. In particular, when c(x, y) = —x ■ y, the convex 
function v provided by Theorem 12.21 solves the classical Monge- Ampere equation 

det(D 2 v(x)) = a.e. 
g{Vv{x)) 

3. The localization argument and proof of the results 

The goal of this section is to prove Theorems 11.31 and 11.41 by showing that the assumptions of 
Theorems 14.31 and 15.31 below are satisfied near almost every point. 

The rough idea is the following: if x is a point where the semiconvex function u is twice differ- 
entiable, then around that point u looks like a parabola. In addition, by looking close enough to 
x, the cost function c will be very close to the linear one and the densities will be almost constant 
there. Hence we can apply Theorem l4.3l to deduce that u is of class C 1 '^ in neighborhood of x (resp. 
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u is of class C k+2 ' a by Theorem 15 .31 if c £ C^ 2 '" and f,g£ Cw*)) wn ich implies in particular that 
T u is of class C 0,/3 in neighborhood of x (resp. T u is of class C k+1 ' a by Theorem 15,31 if c € C^ 2,a 
and f,g£ Cw*)- Being our assumptions completely symmetric in x and y, we can apply the same 
argument to the optimal map T* sending g onto /. Since T* = (Tu)" 1 (see the discussion below), 
it follows that T u is a global homeomorphism of class (resp. T u is a global diffeomorphism of 
class Cj^" 1 '") outside a closed set of measure zero. 
We now give a detailed proof. 

Proof of Theorem 11.31 Let us introduce 
defined as 

u c {y) := 

Then u c is c*-convex, where 

(3.1) c* (y , x) := c(x , y) , and x £ d c *u c (y) <^ y £ d c u(x) 

(see for instance [371 Chapter 5]). 

Being our assumptions completely symmetric in x and y, c* satisfies the same assumptions as c. 
In particular, by Theorem 12. 1| there exists an optimal map T* (with respect to c*) sending g onto 
/. In addition, it is well-known that T* is actually equal to 

T u c(y) = c*-e Wy (Vu c (y)), 

and that T u and T u c are inverse to each other, that is 

(3.2) T u c(T u (x)) = x, T u (T u c(y)) = y for a.e. x £ X , y £ Y 

(see, for instance, Q3 Remark 6.2.11]). 

Since semiconvex functions are twice differentiable a.e., there exist sets X\ C X, Y\ C Y of full 
measure such that (|3.2p holds for every x £ X\ and y £ Y\, and in addition u is twice differentiable 
for every x £ X\ and u c is twice differentiable for every y £ Y\. Let us define 

x' -Xxn^)- 1 ^)- 

Using that T u transports / on g and that the two densities are bounded away from zero and infinity, 
we see that X' is of full measure in X. 

We fix a point x £ X' . Since u is differentiable at x (being twice differentiable), it follows by 
(12. 6p that the set d c u(x) is a singleton, namely d c u(x) = {c-exp s (Vti(x))}. Set y := c-exp i (Vu(x)). 
Since y £Y\ (by definition of X'), u c is twice differentiable at y and x = T u c (y). Up to a translation 
in the system of coordinates (both in x and y) we can assume that both x and y coincide with the 
origin 0. 

Let us define 

u(z) := u(z) - u(0) + c(z, 0) - c(0, 0), 

c(z, w) := c(z, w) — c(z, 0) — c(0, w) + c(0, 0), 

u 5 (w) := u c (w) - u(0) + c(0, - c(0, 0). 

Then n is a c-convex function, u c is its c-conjugate, T a = T u , and T u c = T u c, so in particular 
(Tu)u f = g and (T a c)^ g = f. In addition, because by assumption £ X' , u is twice differentiable 



the "c-conjugate" of u, that is, the function u c : Y 1R 
: sup { - c(x,y) - u(x)}. 



6 G. DE PHILIPPIS AND A. FIGALLI 

at and u c is twice differentiable at = T„(0). Let us define P := D 2 u(0), and M := D xy c(0, 0). 
Then, since c(-, 0) = c(0, •) = and c 6 C 2 , a Taylor expansion gives 

u(z) = ^Pz ■ z + o(|z| 2 ), c(z,w) = Mz ■ w + o(\z\ 2 + \w\ 2 ), 

Let us observe that, since by assumption / and g are bounded away from zero and infinity, by 
(C3) and f|2.10|) applied to u and c we get that det(P), det(M) / 0. In addition (f!T9j) implies that 
P is a positive definite symmetric matrix. Hence, we can perform a second change of coordinates: 
z i — y z : = P 1 / 2 ^, w \— > w := —P~ 1 l 2 M*w (M* being the transpose of M), so that, in the new 
variables, 

(3.3) u(z) := u(z) = ^l^l 2 + °(l^| 2 )> c(z,w) := c(z,w) = —z ■ w + o(\z\ 2 + \w\ 2 ). 
By an easy computation it follows that (!&)«/ = g, where 

(3.4) f{z) := det(P^/ 2 ) f(P~^ 2 z), g(w) := det ((M*)~ 1 P 1 / 2 ) g^M*)' 1 P^ 2 w). 
Notice that 

(3.5) D 22 c(0,0) = £^6(0,0) = n xn, -£>**£((), 0) = Id, D 2 u(0)=ld, 
so, using (|2.10p . we deduce that 

/(0) _ det (D 2 u(0) + £>ggc(0,0)) = x 

3(0) |det(Z>ftBc(0,0))| 

To ensure that we can apply Theorems 14.31 and 15.31 we now perform the following dilation: for 
p > we define 

1 _ 1 

u p {z) := -zu(pz), c p (z,w) := -zc(pz,pw). 
pi pi 

We claim that, provided p is sufficiently small, u p and c p satisfy the assumptions of Theorems 14.31 
and 15.31 

Indeed, it is immediate to check that u p is a c^-convex function. Also, by the same argument 
as above, from the relation (T^)jj/ = g we deduce that T u sends f{pz) onto g{pw). In addition, 
since we can freely multiply both densities by a same constant, it actually follows from (|3.6p that 
( T u p )$f P = 9 P , where 

:= — , o.to := — 

JpV j /(0) Pl ' 3(0) 

In particular, since / and g are continuous, we get 

(3.7) \f p - 1| + \g p - 1| inside £ 3 

as /) — )• 0. Also, by (|3.3|) we get that, for any z,w E B3, 



(3.8) u p (z) = ^\z\ 2 + o(l), c p (5, u;) = -z ■ w + o(l), 



1 An easy way to check this is to observe that the measures fj, := f(x)dx and ^ := g(y)dy are independent of the 
choice of coordinates, hence (|3.4[) follows from the identities 

f(x)dx = f(x)dx, g{y)dy = 
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where o(l) — > as p — > 0. In particular, (|4.9|) and (|4.10|) hold with any positive constants Sq,t]q 
provided p is small enough. 

Furthermore, by the second order differentiability of u at it follows that the multivalued map 
z i — y d~u(z) is differentiable at (see [67\ Theorem 14.25]) with gradient equal to the identity 
matrix (see (I3.3P ). hence 

d~u p (z) C B 7p (z) VzGB 2 , 

where 7 P — > as p — > 0. Since d Cp u p C c p -exp (d~u p ) (by (|2.6p ) and || c p -exp — Id ||oo = o(l) (by 
(PLED), we get 

(3.9) d Cp u p (z) C B Sp (z) Vz£B 3 , 

with 5p = o(l) as p — t- 0. Moreover, the Cp-conjugate of u p is easily seen to be 

u p p (w) = u^{p{M*Y l P l l 2 ui). 

Since u c is twice differentiable at 0, so is u P p . In addition, an easy computation shows that 
D 2 u p p (0) = Id. Hence, arguing as above we obtain that 

(3.10) d Cp u p p (w) cB 5p (w) VweB 3 , 

with 5'p = o(l) as p — > 0. 
We now define 

d:=B u C 2 :=d Cp u p {C{). 

Observe that both C\ and C 2 are closed (since the osubdifferential of a compact set is closed). Also, 
thanks to (|3.9p , by choosing p small enough we can ensure that Bi / 3 C C 2 C B% . Finally, it follows 
from (pUI) that 

(T u )~ 1 {C 2 ) \ C\ C (T^)" 1 ({points of non-differentiability of u p p }), 
and since this latter set has measure zero, a simple computation shows that 

(^ P ) tt (/pld) = 9plc 2 - 

Thus, thanks to (14. 8ft . we get that for any /3 < 1 the assumptions of Theorem 14.31 are satisfied, 
provided we choose p sufficiently small. Moreover, if in addition c G C^. 2 ' a (X x Y), / e C^(X), 
and 5 G Cit'"(^)' t nen a ls° the assumptions of Theorem 15.31 are satisfied. 

Hence, by applying Theorem 14.31 (resp. Theorem 15 .3 j) we deduce that u p G C 1,/3 (l? 1 / 7 ) (resp. 
u p G C k+2,a (B 1 / 9 )), so going back to the original variables we get the existence of a neighborhood Ux 
of x such that u G C^ipk) (resp. u G C k+2 > a (U s )). This implies in particular that T u G C ^^) 
(resp. T u G C k+l ' a {Ux))- Moreover, it follows by Corollary |4.6l that T u {Ux) contains a neighborhood 
of y. 

We now observe that, by symmetry, we can also apply Theorem 14.31 (resp. Theorem 15. 3p to u p . 
Hence, there exists a neig hborhood Vy of y such that T u c G C°' /3 (V S ). Since T u and T u c are inverse 



For instance, this follows by differentiating both relations 

D;c p (5,T Up (5)) = -Vit p (S) and D a c p (T u ^ (w),w) = -\7u p "(w) 

at 0, and using then ([33) and the fact that VT^p (0) = [VT Up (0)]~ 1 and D 2 u p (0) = Id. 
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to each other (see (|3.2p ) we deduce that, possibly reducing the size of U x , T u is a homeomorphism 
(resp. diffeomorphism) between U x and T U {U X ). Let us consider the open sets 

X" := (J U s , Y" := (J T U (Z4), 

and define the (relatively) closed S x := X \ X" , £y := Y \ Y". Since X" D X' , X" is a set of full 
measure, so = 0. In addition, since £y =Y \ Y" C.Y\ T U {X') and T U (X') has full measure 
in Y, we also get that |£y| = 0. 

Finally, since T u : X \ T,x —>Y \ £y is a local homeomorphism (resp. diffeomorphism), by (|3.2I) 
it follows that T u : X \ T,x — > Y \ Xy is a global homeomorphism (resp. diffeomorphism), which 
concludes the proof. □ 

Proof of Theorem \1.4\ The only difference with respect to the situation in Theorem 11.31 is that now 
the cost function c = d 2 /2 is not smooth on the whole M x M. However, even if d 2 /2 is not 
everywhere smooth and M is not necessarily compact, it is still true that the c-convex function 
u provided by Theorem 12.31 is locally semiconvex (i.e., it is locally semiconvex when seen in any 
chart) [131 02] • ^ n addition, as shown in [9j Proposition 4.1] (see also [HI Section 3]), if u is 
twice differentiate at x, then the point T u {x) is not in the cut-locus of x. Since the cut-locus is 
closed and d 2 /2 is smooth outside the cut-locus, we deduce the existence of a set X of full measure 
such that, if xq £ X, then: (1) u is twice differentiable at xq; (2) there exists a neighborhood 
Uxo x Vt u {xq) C M X M of (x ,T u (x Q )) such that c G C°°(U Xo x V Tu(xo) ). Hence, by taking 
a local chart around (xq,T u (xq)), the same proof as the one of Theorem 11.31 shows that T u is 
a local homeomorphism (resp. diffeomorphism) around almost every point. Using as before that 
T u : M — > M is invertible a.e., it follows that T u is a global homeomorphism (resp. diffeomorphism) 
outside a closed singular set of measure zero. We leave the details to the interested reader. □ 

4. C 1 '^ REGULARITY AND STRICT C-CONVEXITY 

In this and the next section we prove that, if in some open set a c-convex function u is sufficiently 
close to a parabola and the cost function is close to the linear one, then u is smooth in some smaller 
set. 

The idea of the proof (which is reminiscent of the argument introduced by Caffarelli in [6] to 
show W 2 ' p and C 2 ' a estimates for the classical Monge- Ampere equation, though several additional 
complications arise in our case) is the following: since the cost function is close to the linear one and 
both densities are almost constant, u is close to a convex function v solving an optimal transport 
problem with linear cost and constant densities (Lemma 14. ip . In addition, since u is close to a 
parabola, so is v. Hence, by [18] and Caffarelli's regularity theory, v is smooth, and we can use this 
information to deduce that u is even closer to a second parabola (given by the second order Taylor 
expansion of v at the origin) inside a small neighborhood around of origin. By rescaling back this 
neighborhood at scale 1 and iterating this construction, we obtain that u is C 1 ^ at the origin for 
some (3 £ (0, 1). Since this argument can be applied at every point in a neighborhood of the origin, 
we deduce that u is C 1 '^ there, see Theorem 14.31 (A similar strategy has also been used in [7] to 
show regularity optimal transport maps for the cost \x — y\ p , either when p is close to 2 or when X 
and Y are sufficiently far from each other.) 

Once this result is proved, we know that d~u is a singleton at every point, so it follows from 
(USD that 

d c u(x) = c-exp x (d~u(x)), 
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see Remark 14,41 below. (The above identity is exactly what in general may fail for general c-convex 
functions, unless the MTW condition holds [31].) Thanks to this fact, we obtain that u enjoys a 
comparison principle (Proposition I5.2| ). and this allows us to use a second approximation argument 
with solutions of the classical Monge- Ampere equation (in the spirit of [6j [27] ) to conclude that u 
is C 2,u in a smaller neighborhood, for some a' > 0. Then higher regularity follows from standard 
elliptic estimates, see Theorem 15.31 

Lemma 4.1. Let C\ and C2 be two closed sets such that 
(4.1) B 1/K CC U C 2 C B K 

for some K > 1, f and g two densities supported respectively in C\ and C2, and u : C\ — > R a 
c-convex function such that d c u(C\) C Bk and (T u )$f = g. Let p > be such that \C\\ = I/0C2I 
(where pC2 denotes the dilation of C2 with respect to the origin), and let v be a convex function 
such that Vfjjlci = lpC2 an d V W = u(0). Then there exists an increasing function uj : R + — > R + , 
depending only K , and satisfying ui(6) > 5 and uj(0 + ) = 0, such that, if 

(4-2) llZ-leJoo + Hs-lcJoo^ 

and 

(4-3) \\c(x,y) + x-y\\ C 2(B K xB K ) < S, 

then 

\W ~ v\\ c o(b 1/k ) < w(<5). 

Proof. Assume the lemma is false. Then there exists Eq > 0, a sequence of closed sets C\, C\ 
satisfying (I4.ip . functions fh, gh satisfying (I4.2p with 5 = 1/h, and costs Ch converging in C 2 to 
—x ■ y, such that 

u h (0) = v h (0) = and \\u h - v h \\ O 0(B 1/K ) > ^o, 
where Uh and Vh are as in the statement. First, we extend Uh an Vh to Bx as 

u h (x) := sup {u h (z)-c h (x,y) + c h (z,y)}, v h (x) := sup {v h (z)+p-(x- z)} . 

zeC^,yed Ch u h {z) zec^,ped-v h {z) 

Notice that, since by assumption d Ch Uh{Ci) C Bk, we have d Ch Uh(Bx) C Bk- Also, (T Uh )$fh = gh 
gives that f fh = f gh, so it follows from (|4.2|) that 

Ph= (\C*\/\C%\y /n ^ 1 ash^oo, 

which implies that d~Vh{Bx) C B Ph x C B2K for h large. Thus, since the C^-norm of Ch is uniformly 
bounded, we deduce that both Uh and Vh are uniformly Lipschitz. Recalling that Uh(0) = %(0) = 
0, we get that, up to a subsequence, Uh and Vh uniformly converge inside Bk to and Voo 
respectively, where 

(4.4) noo (0) = Voo(0) = and \\uoo - 'Uoo||c (s 1/Jr ) > ^0- 

In addition fh (resp. gh) weak-* converge in L°° to some density /oo (resp. goo) supported in Bk. 
Also, since ph — > 1, using f|4.2j) we get that l C h (resp. l ph c^) w eak-* converges in L°° to f^ (resp. 
goo)- Finally we remark that, because of (|4.2p and the fact that C\ D Bi/ K , we also have 

/cxd > lfl 1/Jr - 
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In order to get a contradiction we have to show that u^, = in Bij K . To see this, we apply 
\2>7\ Theorem 5.20] to deduce that both Vmqo and Vt>oo are optimal transport maps for the linear 
cost —x • y sending onto g^. By uniqueness of the optimal map (see Theorem 12. 2\\ we deduce 
that Vf oo = Vmqo almost everywhere inside B\ C spt , hence Mqo = Voc in B\ (since 
Moo(0) = Woo(0) = 0), contradicting fljgp . □ 

Here and in the sequel, we use Af r (E) to denote the r-neighborhood of a set E. 

Lemma 4.2. Let u and v be, respectively, c-convex and convex, let D G M nx ™ be a symmetric 
matrix satisfying 

(4.5) Id/K < D < Kid 

for some K > 1, and define the ellipsoid 

E(xq, h) := {x : D{x — xq) • (x — xq) < h], h > 0. 
Assume that there exist small positive constants e, 5 such that 

( 4 - 6 ) \\ v ~ u \\c°(E(x ,h)) <£, \\c + X- y\\c2(E(x ,h)xd c u(E(x ,h)) < 5 - 

Then 

(4.7) d c u(E(xo,h-y/i)) cN K , {S+vre) (dv(E(x ,h))) V0<e<h 2 <l, 

where K' depends only on K . 

Proof. Up to a change of coordinates we can assume that xq = 0, and to simplify notation we set 
Eh '■= E(xQ,h). Let us define 

v(x) := v(x) + e + 2yfe(Dx ■ x — h), 

so that v > u outside Eh, and v < u inside E h _^. Then, taking a c-support to u in E h _^ (i.e., a 
function C XjV as in (|2.3p . with x G E h _^ and y G d c u(x)), moving it down and then lifting it up 
until it touches v from below, we see that it has to touch the graph of v at some point x G Eh- in 
other words 

d c u{E h _^) C d c v{E h ). 

By (|4.5p we see that diam^, < y/Kh, so by a simple computation (using again (|4.5|) ) we get 

Thus, since d c v(Eh) C c-exp(3 _ -y(£'/ l )) (by (|2.6|) ) and || c-exp — Id ||c;o < 5 (by (|4.6|) ). we easily 
deduce that 

dcu{E h _^) C M RI{s+Vhl) {d-v{E h )), 
proving (|4.7p . □ 
Theorem 4.3. Xei Ci and C2 6e two closed sets satisfying 

B 1 / 3 C Ci,C 2 C B 3 , 



Even if v is not c-convex, it still makes sense to consider his e-subdifferential (notice that the c-subdifferential 
of v may be empty at some points). In particular, the inclusion d c v(x) C c-exp x (d~ v(x)) still holds. 
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let f,g be two densities supported in C± and C2 respectively, and let u : C\ — > R be a c-convex 
function such that d c u{C\) E B3 and (T u )$f = g. Then, for every (3 E (0, 1) there exist constants 
Sq,7]q > such that the following holds: if 



(4.8) 

(4.9) 

and 

(4.10) 



3-Cilloo + ||<7 - lCalloo < 



\c(x,y) + x- y\\ C 2 



{B 3 xB 3 ; 



<S , 



u 



< 



C°(B 3 ) 



v \x\ 

T 1 



then u E C 1 '^(B 1/7 ). 

Proof. We divide the proof into several steps. 

• Step 1: u is close to a strictly convex solution of the Monge Ampere equation. Let v : R n — > R 
be a convex function such that Vujjld = lpc 2 with p = (|Ci|/|C2|) 1//n (see Theorem 12 .2p . Up to a 
adding a constant to v, without loss of generality we can assume that v(0) = u(0). Hence, we can 
apply Lemma 14. II to obtain 

( 4 - n ) \\ v ~ u \\c°(B 1/3 ) < "(fa) 

for some (universal) modulus of continuity uj : R + — > R + , which combined with (|4,10p gives 

< r/o + u(S ). 

C°(B 1/3 ) 

Also, since fc 1 f = Jc 9i ^ follows easily from (|4.8p that \p — 1| < 3^o- By these two facts we 
get that d~v(Biu) C Bj/24 C pC2 provided 5q and r/o are small enough (recall that v is convex 
and that B 1 / 3 C C2), so we can apply [181 Proposition 3.4] to deduce that v is a strictly convex 
Alexandrov solution to the Monge- Ampere equation 

(4.12) detD 2 v = l in B 1/A . 

In addition, by a simple compactness argument, we see that the modulus of strict convexity of v 
inside B 1 u is universal. So, by classical Pogorelov and Schauder estimates, we obtain the existence 
of a universal constant Kq > 1 such that 

(4.13) \\v\\ cHBi/5 ) < K , U/K <D 2 v<K ld in B 1/5 . 
In particular, there exists a universal value h > such that, for all x E Buy, 

Q(x,v, h) := {z : v(z) < v(x) + Vtj(x) • (z - x) + h} EE B 1 / 6 Vh<h. 

• Step 2: Sections of u are close to sections of v. Given x E B 1 /> 7 and y E d c u(x), we define 

S(x,y,u,h) := {z : u(z) < u(x) - c(z,y) + c(x,y) + h}. 
We claim that, if 5q is small enough, then for all x E -B1/7, y E d c u(x), and h < h/2, it holds 

(4.14) Q(x,v,h- K iy /u3(6 )) E S(x,y,u,h) E Q(x, v, h + K iy /u(So)), 
where K\ > is a universal constant. 
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Let us show the first inclusion. For this, take x £ -B1/7, y G d c u(x), and define 

Px ■= -D x c(x,y) £ d~u(x). 

Since v has universal C 2 bounds (see (|4.13|) ) and u is semi-convex (with a universal bound), a 
simple interpolation argument gives 

(4.15) \p x - Vv(x)\ < K'yJWu - v\\ c o {Bi/5) < K'^J(5o) VxeB 1/T . 
In addition, by (fill . 

(4.16) \y~Px\ < \\D x c + Id\\ c o (b 3 xB 3 ) < 5 o, 
hence 

(4.17) + c(z,y)| < |z-p x - « ■ y\ + |z- 2/ + c(z,y)\ < 28 \/x,zeB 1/7 . 
Thus, if z € Q(x, v,h- Kiy/u(6 )), by gZTTJ) , pi5j> . and (|4TT7D we get 



u(z) < v(z) + w(5o) < v(x) + Vv(x) ■ (z - x) + h - K\y/u(8o) + uj(8 ) 

< u(x) + p x -z-p x -x + h- Ki\Juj(5q) + 2u){8q) + 2K' \/ uj(5q) 

< u(x) - c(z, y) + c(x, y) + h- K x a/w(<5 ) + 2w(<5 ) + 2KV w (^o) + 45 

< u(x) - c(z, y) + c(x, y) + h, 

provided K\ > is sufficiently large. This proves the first inclusion, and the second is analogous. 

• Step 3: Both the sections of u and their images are close to ellipsoids with controlled eccentricity, 
and u is close to a smooth function near xq. We claim that there exists a universal constant 
K2 > 1 such that the following holds: For every 770 > small, there exist small positive constants 
ho = ho(r]o) and So = 5o(ho,r]o) such that, for all xo £ -B1/7, there is a symmetric matrix A satisfying 

(4.18) Id /K 2 < A < K 2 Id, det(A) = 1, 

and such that, for all yo 6 d c u(xo), 



A IB 



(4.19) 

Moreover 
(4.20) 



^j^(xo)) C S(xo,yo,u,h ) 



C A B 



/8W 



Oo) 



A 1 OVvb^) c d c u(S(xo,y ,u,ho)) c A 1 (#^(1/0)) • 



u-a 



-\A'\x-xo)\ 



C° A B 



^0 



(*o))) 



< mho, 



where C Xoyo is a c-support function for u at xo, see (|2.3p . 

In order to prove the claim, take ho <C h small (to be fixed) and 80 <C ho such that Kiy/u)(8a) < 
ho/2, where K\ is as in Step 2, so that 

(4.21) Q(x ,v, h /2) C S(x , y , u, h ) C Q(x , v, 3/i /2) CC B 1/6 . 
By (|4.13p and Taylor formula we get 

(4.22) v(x) = v(x ) + Vv(x ) ■ (x - x ) + -D 2 v(x )(x - x ) ■ (x - x ) + 0(\x - x \ 3 ), 
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so that defining 

(4.23) E(x ,h ) := lx : ^D 2 v(x )(x - x ) • (x - x ) < /i | 

and using (|4.13p . we deduce that, for every ho universally small, 

(4.24) E(x , ho/2) C Q(x , v, h ) c E(x , 2h ). 
Moreover, always for ho small, thanks to (|4.22p and the uniform convexity of v 

(4.25) Vv(E(x , h )) C E*(Vv(x ), 2h ) C Vv(E(x , 3h )) 
where we have set 

E*(yM) ■= |y : \[D 2 v(y)Y\y-y)-{y-y) </i j. 

By Lemma I3~2j (I4.24p . and (|4.25p applied with 3/io in place of ho, we deduce that for 5o <S ho <S 

(4.26) d c u(S{x ,yo,u,h )) cM K „ y ^(Vv(E(x ,ah ))) C E*(Vv(x ), 7h ). 

Moreover, by (|4.15p . if yo E d c u(xo) and we set p XQ := —D x c{xo,yo), then 

\y - Vv(x Q )\ < \p XQ - Vv(x )\ + \\D x c + ld\\ c o^ B3XB ^ < \/u)(6 ) + S - 
Thus, choosing 5o sufficiently small, it holds 

(4.27) E*(Vv(x Q ),7h )cE*(y ,8h ) Vy G d c u(x ). 
We now want to show that 

E*(yo,h /8) C d c u(S(x ,yo,u, h )) Vy G d c u(x ). 
Observe that, arguing as above, we get 

(4.28) E*(y , ho/ 8) C E*(Vv(x ), ho/7) Vy G 9 c n(x ) 
provided Jo is small enough, so it is enough to prove that 

E*(Vv(x ),h /7) c d c u(S(x ,yo,u,h )). 
For this, let us define the c*-convex function u c : B% — > R and the convex function t> * : B3 — > R 



(2/):= sup { - c(x,y) - u(x)}, v*(y):= sup {a; • y - v(x)} 

xeB 1/5 ieBi/5 



(see (|3.ip ). Then it is immediate to check that 

(4.29) \u c -v*\< lo(5o) +5 < 2u(5 ) on B 3 . 

Also, in view of (I4.13p . v* is a uniformly convex function of class C 3 on the open set Vf (-B1/5) 
addition, since 

(4.30) F C d c u(d c *u c (F)) for any set F, 
thanks to (|4.2ip and (|4.24p it is enough to show 

(4.31) d c *u c (E*(Vv{x ),h /7)) C E(x ,ho/A). 
For this, we apply Lemma 14.21 to u c and v* to infer 

d c *u c (E^Vv(x ), ho/7)) C Af Kf „s^{Vv*(E*(Vv(xo), ho/7))) C E(x ,h /4), 
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where we used that 



\7v* = [Vv]" 1 and D 2 v*(Vv(x )) = [D 2 v (xq)}' 1 . 



Thus, recalling ()4.26p . we have proved that there exist ho universally small, and 6$ small depending 
on ho, such that 



(4.32) 



E*(Vv{x ),h /7) C d c u{S{x ,yo,u,h )) C E*(Vv(x ),7h ) Vx G B 1/7 . 



Using g20, (gUZD , and ([4T2"g|) . this proves (|4TT9|) with 4 := [D 2 v(x )}' 1/2 ■ Also, thanks 

to (j4T2D and (|4T3]), (g38j) holds. 

In order to prove the second part of the claim, we exploit ([4.11 p . (|4.9p . (|4.16|) . (|4.15|) . (|4.22|) . and 
(ETTKI) (recall that C Xom is defined in JZI]) and that A = [D 2 v(x )}~ 1/2 ): 



u-C, 



XQ,V0 ( X 



1 



C°(E(x ,8h )) 



u - C Xo , yo - v(x )(x - x Q ) ■ (x - x ) 

1 C°(E(x ,8h )) 
< 2\\u - u||cfO(B( a . ,8fco)) + M x ,yo) + X ■ yo\\ C °(E(xo,8h )) + \\ c ( x O,Vo) + ^0 • Vo\\c°{E(xo 



8ho)) 



+ 



+ 



\{VO-Pxo) ■ ( x ~ x o)\\ c o {E{ 



xo,8ho)) 



+ 



( Pxo -Vv(x )) ■ ( x ~ x o)\\ c o (E{x ^ 



Sho)) 



1 



v - v(x ) - Vv(x ) ■ (x - x ) - -D v(x )(x - x ) ■ (x - Xq) 



< 2u(5 ) + 35 + K'^/uj(5o) + K i^V^o <Voh , 



C0(E(x ,8h )) 



where the last inequality follows by choosing first ho sufficiently small, and then 5o much smaller 
than ho- 
rn Step 4-' A first change of variables. Fix xo £ Bxn, yo G d c u(xo), define M := —D xy c(xo,yo), and 
consider the change of variables 

{x := x — xq 
y:=M- x {y-yo). 

Notice that, by ([4~9|h it follows that 

(4.33) |M-Id| + |M -1 -Id| <35 
for 5o sufficiently small. We also define 

c(x, y) := c(x, y) - c(x, yo) - c(x , y) + c(x , yo), 
u(x) := u(x) - u(x ) + c(x, y ) - c(x ,y ), 
u c (y) := u c (y) - u c (y ) + c(x ,y) - c(x , yo). 
Then u is c-convex, u c is c*-convex (where c*(y,x) = c(x,y)), and 

(4.34) c(-,0) = c(0,-) = 0, D xg c(0,0) = -Id. 
We also notice that 

(4.35) 8- c u{x) = M" 1 (d c u{x + xo) - y ) . 
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Thus, recalling (|4.19p . and using (|4.33|) and (|4.35|) . for <5o sufficiently small we obtain 
(4.36) 



A B 



A' 1 IB 



VW9 
y/h /9 



C 5(0, 0, u, h ) CA[B 



c M~ A 



i a-i 



B 



1-1 



B 



C d 5 u(S(0,0,u,h )) C M~ 1 A~ 1 [B^j C A- 
Since (T M )jj/ = g, it follows that T s = c-exp(Vu) satisfies 

(2fl)|j/ = 5, with /» := f(x + x ), := det(M) g(My + y ) 

(see for instance the footnote in the proof of Theorem 1 1.3|) . Notice that, since \M — Id | < 5q (by 
(|4U1) L we have | det(M) - 1| < (1 + 2ra)<5 (for 5 small), so by flMD we get 

(4-37) 11/ - l Cl -x ||oo + \\9 ~ lM-i(Ca- M )lloo < 2(1 + n)S . 

• Step 5: A second change of variables and the iteration argument. We now perform a second 
change of variable: we set 



(4.38) 
and define 



nf • \ 1 rv* 

Jj . 7== -Tl Jb 

V «0 

V := VGA*' 



ci(x,y) := — c(-\fhoAx, \fh~oA 

u\(x) := — u( \fho~Ax), 
ho 

u?(y) ,= Lu\^hoA~ l y). 
ho 



We also define 



fi(x) := KVhoAx), g x {y) := g^^A^y). 
Since det(^4) = 1 (see (]4.18p ). it is easy to check that (T Ul )$fi = g\ (see the footnote in the proof 
of Theorem [L3]). Also, since (\\A\\ + ^/ho~ < 1, it follows from (147371) that 

(4.39) \ft-l\ + \gi - 1| < 2(1 + n)5 inside B 3 . 
Moreover, defining 

C{ 1] := 5(0,0,^,1), eg* :=d clUl (5(0, 0,«i,l)), 
both and are closed, and thanks to (14.36H 

(4.40) ^cClVfcBa. 

Also, since (T Ul )$fi = g\, arguing as in the proof of Theorem 11.31 we get 

(r ui )tt(/il c (i)) = (9i1 c (d), 

and by (Q9l) 

||/il (i) - l (i) ||oo + ||5il r (i) - l r (i)||oo < 2(1 + n)5 . 

u l u l u 2 u 2 

Finally, by (|4.34p and (|4.20p , it is easy to check that 

1, 



\ci{x,y) + x •y||c2(B 3 xB 3 ) ^ <*o, 



ui \x\ 2 

2 I I 



< 



C°(B 3 ) 
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This shows that U\ satisfies the same assumptions as u with 5q replaced by 2(1 + n)6o. Hence, up 
to take 5q slightly smaller, we can apply Step 3 to ui, and we find a symmetric matrix A± satisfying 

Id /K 2 <At<K 2 Id, det(Ai) = 1, 



Ai B 



/ho/8 



C 5(0,0,«i, /t ) C Ax IB 



/Mo > 



V (s^J C 9 cl m (5(0, 0, Ul ,/io)) C V (5 

1| „-l-|2 

ui - -\A 1 x\ 



/8ho 



< mho- 

C°(A 1 (B(0,V8h ~)) 
(Here K 2 and ho are as in Step 3.) 

This allows us to apply to u\ the very same construction as the one used above to define u\ from 
u: we set 

C2{x,y) := —ciUfhoAix, y/hoA^y), u 2 {x) := —u 1 (y/h ~Aix), 
no no 

so that (T U2 )jj/ 2 = g 2 with 

f2{x) := h{\fh^Aix), g 2 {y) ■= giy/hoA^y). 

Arguing as before, it is easy to check that u 2 , C2, f2,92 satisfy the same assumptions as u\, c±, fi,gi 
with exactly the same constants. 

So we can keep iterating this construction, defining for any k £ N 

Ck+i(x,y) := —c k (\/%A k x,y/h^A7 1 y), u k+x (x) := -^-u k (y/h^A k x), 
ho no 

where A k is the matrix constructed in the A:-th iteration. In this way, if we set 

M k :=A k -...-A 1 , Vfc>l, 
we obtain a sequence of symmetric matrices satisfying 
(4.41 ) Id /Kl < M k < Kl Id, det (M fc ) = 1 , 

and such that 



(4.42) 



M ^ ( S (ho/8)*/a) C 5(0,0,^,^) C M k (B (8ho)k/ ^ . 



• S^ep 6: C 1 ^ regularity. We now show that, for any (3 £ (0, 1), we can choose ho and <5o = ^o(^-o) 
small enough so that u\ is C 1 '^ at the origin (here u\ is the function constructed in the previous 
step). This will imply that u is C 1 '^ at xo with universal bounds, which by the arbitrariness of 
xo € B 1/7 gives u £ C 1 ^(B l/1 ). 

Fix P G (0, 1). Then by (IOTP and (jl^2l) we get 



B, s.k C S(0,0,ui,hn) C B, o 



(4.43) 

so defining ro := y/h~o/(y/8K2) we obtain 

||«i|bo(s ») < = (v^^ro) 2 " < r Q 1+/3)fc 



provided /io (and so ro) is sufficiently small. This implies the C 1 '^ regularity of u\ at 0, concluding 
the proof. □ 
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Remark 4.4 (Local to global principle). If u is differentiable at x and c satisfies (CO)-(Cl), then 
every "local c-support" at x it is also a "global support" at x, that is, d c u(x) = c-exp x (d~ u(x)) . 
To see this, just notice that 

7^ d c u(x) C c-exp x (d~ u(x)) = {c-exp x (Vu(x))} 

(recall (j2.6|) ) . so necessarily the two sets have to coincide. 

Corollary 4.5. Let u be as in Theorem \4-3\ Then u is strictly c-convex in Bij 7 . More precisely, 
for every 7 > 2 there exist 770, <5o > depending only on 7 such that, if the hypothesis of Theorem 
\4-3\ are satisfied, then, for all xo G -B1/7, yo G d c u(xo), and C XOtVO as in (12.3p . we have 

(4.44) inf {u - C x ^ yo ] > cor 7 Vr< dist(x , 8B 1/7 ), 

8B r (x ) 

with Co > universal. 

Proof. With the same notation as in the proof of Theorem I4.3| it is enough to show that 

inf u\ > r 1///3 , 

8B r 

where u\ is the function constructed in Step 5 of the proof of Theorem l4.31 Defining ^0 := K2y/8ho, 
it follows from (|4.43|) that 

inf Ul > h k = ( Q o/(V8K 2 )) 2k > of, 

provided ho is small enough. □ 

A simple consequence of the above results is the following: 
Corollary 4.6. Let u be as in Theorem \4.3\ then T u {B^p) is open. 

Proof. Since u S C 1 '^(S 1 / 7 ) we have that T u (B 1 / 7 ) = d c u(Bi/ 7 ) (see Remark l4.4p . We claim that 
it is enough to show that if yo G d c ii(B 1 / 7 ), then there exists e = e(yo) > small such that, for all 
\y — 2/0 1 < £, the function u(-) + c(-, y) has a local minimum at some point x G -B1/7. Indeed, if this 
is the case, then 

Vu(x) = -D x c(x,y), 
and so y G d c u(x) (by Remark l4.4p . hence B £ {yo) C T u (B±/ 7 ). 

To prove the above fact, fix r > such that B r (xo) C -B1/7, and pick x a point in i? r (xo) where 
the function u(-) + c(-,y) attains its minimum, i.e., 

x G argmin \u(x) + c(x,y)\. 

Br(x ) 

Since, by (|Q4) . 

min |ii(x) + c(x, y)} > min {u(x) + c(x, yo)} — e||c||c<i 

x£dB r (xo) x£dB r (xo) 

> u(x ) + c(x , yo) + c r 7 - e||c|| c i , 

while 

u(x ) + c(x , y) < c(x ,yo) + it^o) + e\\c\\c^ , 
choosing e < ^tot 1 r7 we obtain that x G B r (xo) C -B1/7. This implies that x is a local minimum 
for u(-) + c(-,y), concluding the proof. □ 
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5. Comparison principle and C 2,a regularity 
We begin this section with a change of variable formula for the c-exponential map. 
Lemma 5.1. Let Q be an open set, v S C 2 (f2), and assume that S7v(Q) C Domc-exp and that 

D 2 v(x) + D xx c[x, c-exp^ 
Then, for every Borel set A C Q, 

|c-exp(Vt;(A))| < f det { D M%) + Dxxcjx, c-e Wx (Vv(x)))) ^ 
J A |det(D^c(x,c-exp a; (Vf(x))))| 

In addition, if the map x i— >■ c-exp x (Vv(x)) is injective, then equality holds. 

Proof. The result follows from a direct application of the Area Formula [12^ Section 3.3.2, Theorem 
1] once one notices that, differentiating the identity 

Vv(x) = — D x c[x, c-exp a ,(Vv(x))) 
(see (|2.5p ). the Jacobian determinant of the C 1 map x h-» c-exp x (Vf (x)) is given precisely by 

det ( y D 2 v{x) + D xx c(x, c-exp a ,(V?;(x)))) 
\det(D xy c(x, c-exp x (Vi)(x)))) | 

□ 

In the next proposition we show a comparison principle between C 1 c-convex functions and 
smooth solutions to the Monge-Ampere equation. ^ As already mentioned at the beginning of 
Section [J] (see also Remark 14. 4p . the C 1 regularity of u is crucial to ensure that the c-subdifferential 
coincides with its local counterpart c-exp(<9~-u). 

Here and in the sequel, we use co[E] to denote the convex hull of a set E. Also, recall that N r {E) 
denotes the r-neighborhood of E. 

Proposition 5.2 (Comparison principle). Let u be a c-convex function of class C 1 inside the set 
S := {u < 1}, and assume that u(0) = 0, E\ij£ C S C Bx, and that X?u(S) CC Domc-exp. Let 
f,g be two densities such that 

(5-1) ll/Ai - l||c°(S) + II5A2 - M\c°(T u (S)) < £ 

for some constants Ai,A2 G (1/2,2) and e E (0,1/4), and assume that (T u )$f = g. Furthermore, 
suppose that 

(5-2) \\c + x-y\\ C 2(B K xB K ) < S. 

Then there exist a universal constant 7 G (0, 1), and 5i = b~i(K) > small, such that the following 
holds: Let v be the solution of 

(det(D 2 v) = A1/A2 in M 5 -,(co[S]), 
\v = l on d(AMco[S])). 

Then 

(5.3) \\u — u||co(5) < Ck \£ + o" 7 /™^ provided 5 < 5±, 

4 A similar result for the case c(x, y) — \x — y\ p appeared in [7J Theorem 6.2]. Here, however, we have to deal with 
some additional difficulties due to the fact that the c-exponential map is not necessarily denned on the whole R™. 
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where Ck is a constant independent of \\, A2, £, and 5 (but which depends on K). 

Proof. First of all we observe that, since u(0) = 0, u = 1 on dS, S C Bk, and ||c + x ■ y\\c 2 (B K ) — 
8 <C 1, it is easy to check that there exists a universal constant a\ > such that 

(5.4) \D x c(x,y)\ > a\ Vx G dS, y = c-exp x (S/u(x)). 

Thanks to (|5.4p and (|5.2p . it follows from the Implicit Function Theorem that, for each x G dS, 
the boundary of the set 

E x := {z £ B K : c(z, y) - c(x, y) + u(x) < l} 

is of class C 2 inside B^, and its second fundamental form is bounded by CrS, where Ck > 
depends only on K. Hence, since S can be written as 

S := p| E x , 

x£dS 

it follows that 

S is a (CK<5)-semiconvex set, 
that is, for any couple of points xo,x\ G S the ball centered at xi/ 2 ■= (xo + x\)/2 of radius 
Ck5\x\ — xq\ 2 intersects S. Since S C Bk, this implies that cofS"] C Nc' K s(S) f° r some positive 
constant C' K depending only on K. Thus, for any 7 G (0, 1) we obtain 

Since v = 1 on d(Ns7 (co[S*])) and A1/A2 G (1/4,4), by standard interior estimates for solution of 
the Monge- Ampere equation with constant right hand side (see for instance [HJ Lemma 1.1]), we 
obtain 

(5.5) osc s v<C K 

(5.6) 1 - C' K ^l n <v<l ondS, 

(5.7) D 2 v > 5"// T Id /C K in co[5], 

for some r > universal, and some constant C K depending only on K. 
Let us define 

v + ■= (1 + 4e + 2y/S)v -As- 2^6, v~ := (1 - 4e - ^6/2)v + 4e + y/5/2 + 2C' K 5 l/n . 

Our goal is to show that we can choose 7 universally small so that v ~ > u > v + on S. Indeed, if 
we can do so, then by (15. 5p this will imply (I5.3P , concluding the proof. 

First of all notice that, thanks to (|5.6p . v~ > u > v + on OS. Let us show first that v + < v. 

Assume by contradiction this is not the case. Then, since u > v + on dS, 

:= {u< v + ) CC S. 

Since v + is convex, taking any supporting plane to v + at x G Z, moving it down and then lifting 
it up until it touches u from below, we deduce that 

(5.8) X7v+(Z) C Vu{Z) 
(recall that both u and v + are of class C 1 ), thus by Remark | 

(5.9) |c-exp(V?; + (Z))| < \T U (Z)\. 
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We show that this is impossible. For this, using (|5.7p and choosing 7 := r/4, for any x 6 Z we 
compute 

D 2 v + (x) + D xx c(x, c-exp x {Vv + (x))) > (1 + + 4e)D 2 v + V6D 2 v -5 Id 

> (1 + V5 + 4e)D 2 v + (5 1/4 /C£ - 5) Id 

> (1 + y/S + 4e) J D 2 t>, 

provided 5 is sufficiently small, the smallness depending only on K. Thus, thanks (|5.2p we have 

det (D 2 v+(x) + D xx c(x,c-exp x (Vv+(x)))) > det ((1 + y/S + 4e)D 2 v) 
\det(D xy c(x, c-exp 3 ,(Vu+(a;)))) | 1 + 5 

> (1 + Vd + 4e) n (l-25)\ 1 /X 2 

> (l + 4ne)Ai/A 2 . 

In addition, thanks (|5,7p and (|5.2p . since 

§1 /t = §1/4 ^ £ 

we see that 

£> 2 t> + > \\D xx c\\ c o( Bl<xBK )Id inside co [5]. 

Hence, for any x, z £ Z , x ^ z and y = c-exp x (Vv + (x)) (notice that c-ex.p x (Vv + (x)) is well-defined 
because of (|5.8p and the assumption Vu(S) CC Domc-exp), it follows 

v + {z) + c(z, y) > v + {x) + c(x, v) + \J a ( d2v+ ( tz + 0- ~ + D xxc{tz + (1 - t)x, y)) [z - x, z - 
> f + (x) + c(x,y), 

where we used that Vv + (x) + D x c(x,y) = 0. This means that the supporting function z 
—c(z, y) + c(x, y) + v + (x) can only touch v + from below at x, which implies that the map 
c-exp x (Vu + (x)) is injective. Thus, by Lemma IBTTI we get 

(5.10) I c-exp(Vu + (Z))| > (1 + 4ne)X 1 /X 2 \Z\. 

On the other hand, since u is C , it follows from (T u )»f = g and (|5.ip that 

J z g{T u (x)) A 2 (l-e) 

This estimate combined with (|5.10|) shows that (|5.9|) is impossible unless Z is empty. This proves 
that v + < u. 

The proof of the inequality v~ < u follows by the same argument expect for a minor modification. 
More precisely, let us assume by contradiction that W := {u > v~} is nonempty. In order to apply 
the previous argument we would need to know that Vv~(W) C Domc-exp. However, since the 
gradient of v can be very large near dS, this may be a problem. 

To circumvent this issue we argue as follows: since W is nonempty, there exists a positive 
constant p, such that u touches v~ + p, from below inside 5. Let E be the contact set, i.e., 
E := {u = v~ + ji}. Since both u and v~ are C , Vw = Vu~ on E. Thus, if 77 > is small 
enough, then the set := {u > v~ + Ji — 77} is nonempty and Vv~(W r) ) is contained in a small 
neighborhood of Vu(W v ), which is compactly contained in Domc-exp. At this point, one argues 
exactly as in the first part of the proof, with W„ in place of Z, to find a contradiction. □ 
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Theorem 5.3. Let u, /, g, r/o, So be as in Theorem \4-^[ and assume in addition that c E C k,a {B^ x 

1/3 



B%) and f,g E C k,a (Bi/ 3 ) for some k > and a G (0, 1). There exist small constants 771 < 770 and 



Si < So such that, if 

(5.H) ||/-l Cl ||oo + |b-lc 2 ||oo <S 1 

(5-12) \\c(x,y) + x- y\\ C 2(B s xB a ) < Si 
and 



(5.13) 



1, ,2 

u~- 2 \x\ 



< m, 

C°(B 3 ) 



then u G C k+2 ' a (B 1/9 ). 

Proof. We divide the proof in two steps. 



Step 1: C > regularity. Fix a point xq G B 1 / 8 , and set yo ■= c-exp Xo (Vu(xo))- Up to replace u 



(resp. c) with the function u\ (resp. ci) constructed in Steps 4 and 5 in the proof of Theorem 14.3 
we can assume that u > 0, u(0) = 0, that 

S h := S(0,0,u,h) = {u< h}, 

and that 

(5.14) D xy c(0,0) = -U. 

Under these assumptions we will show that the sections of u are of "good shape", i.e., 

(5.15) B^CS.CB^ Vh<h u 

for some universal hi and K. Arguing as in Step 6 of Theorem 14.31 this will give that u is C ' l at 
the origin, and thus at every point in Bi/ 8 . 

First of all notice that, thanks to f|5. 13[) . for any hi > we can choose rj\ = f?i(/ii) > small 
enough such that ()5. 15j) holds for with K = 2. Hence, assuming without loss of generality that 
8\ < 1, we see that 

B Vhl/z <= A/^^cofS^]) C B zVK - v 
where 7 is the exponent from Lemma 15.21 Let vi solve the Monge- Ampere equation 

f det(DV) = /(0)/ 5 (0) mM s ^(co[S hl ]), 
\ vi = hi on 9J\f s y^(co[S hl ]). 

Since -B1/3 C N ~p v ^r(co[5'/ ll ]) / '\fh\ C B3, by standard Pogorelov estimates applied to the function 

vi{yfh~lx)/hi (see for instance [261 Theorem 4.2.1]), it follows that |L>V(0)| < M, with M > 
some large universal constant. 

Let /ifc := hi2~ k and define X > 3 to be the largest number such that any solution w of 

(5.16) ( = * >M°) m ^ with B 1/Jf CZC% 
[ w = 1 on oZ, ^ 
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satisfies \D 2 w(0)\ < M + 1. We prove by induction that (|5,15p holds with K = K. 

If h = hi then we already know that (j5.15[) holds with K = 2 (and so with K = K). 

Assume now that (|5.15|) holds with h = h k and K = K, and we want to show that it holds with 
h = hk+x- For this, for any k G N we consider u k the solution of 

/ det(D 2 v k ) = f(0)/g(0) in JV^(co[S h J), 
\ v k = h 1 2~ k ondN qvTrk (co[S hk }), 

where 

5 k := \\c{x,y) + x ■ y\\ C 2 [ShkXTu{Shk)) < S v 

Let us consider the rescaled functions 

u k (x) := u(y/h^x)/h k , v k {x) := v k (^Jh~ k x) /h k . 

Since by the inductive hypothesis B 1 i R C S k := {u k < 1} C B R , we can apply Lemma [5.21 to 
deduce that 

■I'M ^ n-(x. 



(5.17) \\u k -v k \\ c0{ g k) < C R [ osc/+ osc 9 + 8% ) <C R (8i + 8{ 

This implies in particular that, if 5i is sufficiently small, B 1 ^ 2 R) c {^fc — 1} c ^2A- By standard 
estimates on the sections of solutions to the Monge-Ampere equation, the shapes of {v k < 1} and 
{v k < 1/2} are comparable, and in addition sections are well included into each other \26\ Theorem 
3.3.8]: there exists a universal constant L > 1 such that 

B i/{LR) C {vk < 1/2} C B LR , dist({i; fc < 1/4}, 8{v k < 1/2}) > 1/(L#). 

Using again (|5,17p we deduce that, if 8\ is sufficiently small, 

B 1/{2 lk) C {u fc < 1/2} C B 2L #, dist({n fe < l/4},3{u fc < 1/2}) > 1/(2LA) 

so, by scaling back, 

( 5 - 18 ) B ^/{2LK) C S ^ + l C B %LRy/K£L' ^t(S hk+2 , dS hk+1 ) > VT k /(2LK). 

This allows us to apply Lemma 15.21 also to u k +i to get 

(5-19) -v k+1 \\ CO ( S j < C 2lR ( osc /+ osc s + 

We now observe that, by (|5.15|) and the C 1 ^ regularity of u (see Theorem 14. 3|) . it follows that 

diam(5 /l J + diam(T u (S hk )) < Ch^ 2 , 
so by the C 0,a regularity of / and g, and the C 2,a regularity of c, we have (recall that 7 < 1) 

(5.20) osc / + osc g + 5l' n < C'hZ, a := ^ 

s hfe r u (S h .) k 2n 



^The fact that K is well defined (i.e., 3 < K < 00) follows by the following facts: first of all, by definition, M is 
an a-priori bound for |D 2 ui(0)| whenever w is a solution of (|5. 1 6p with B1/3 C Z C -B3, so A" > 3. On the other hand 
K < \fM + 1, since the function 

2 

uS := (M + l)xj + -^j— +xl + ...+x'i 
is a solution of (pHo) with B 1/v ^t+t C Z C B^jrf+T and l-° 2 w(0)| = 2(M + 1). 
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Hence, by (foTTTD and (EnS]) . 

\\uk - Vk\\c°(s k ) + Nfe+i ~ ^fc+i|lco(5 fc+ i) < C (C R + C 2hR ) h%, 
from which we deduce (recall that h k = 2h k+ i) 

\\vk ~ v k+1 \\ c o {Shk+i) < \\v k - u\\ c o {Shk) + ||u - v k+1 \\ c o {Shk+i) 

= h k \\u k — Vk\\c°(s k ) + %+ill^fc+i ~~ Vk+i\\c°(s k+1 ) 
<C{C R + C 2Lk )h\ + °. 

Since v k and v k+ \ are two strictly convex solutions of the Monge Ampere equation with constant 
right hand side inside Sh k+1 , and since Sh k+2 is "well contained" inside Sh k+1 , by classical Pogorelov 
and Schauder estimates we get 

(5-21) \\D 2 v k - D 2 v k+1 \\ c0(Shk+2) < C' R hl 

(5-22) \\D\ k - D 3 v k+1 \\ C o iShk+2) < C' R hl- 1/2 , 

where C'^ is some constant depending only on K. By (|5.2ip applied to Vj for all j = 1, . . . , k (this 
can be done since, by the inductive assumption, (|5.15p holds for h = hj with j = 1,. . . ,k) we 
obtain 

k 

\D 2 v k+1 (0)\ < \D 2 vi(0)\ + \D\{0) - D 2 v j+l (0)\ 

i=i 

k 

<M + C' R hl 2 ~ ja 
C- 

< M + Z—bK < M + 1, 

1 - 2~ CT 1 ~ 

provided we choose h\ small enough (recall that h k = h\2~ k ). By the definition of K it follows 
that also Sh k+1 satisfies (|5.15p . concluding the proof of the inductive step. 

• Step 2: higher regularity. Now that we know that u £ C 1,1 (i?i/g), Equation (|2.10p becomes 

uniformly elliptic. So one may use Evans-Krylov Theorem to obtain that u € C 2 ^ (-B1/9) for some 
a' > 0, and then standard Schauder estimates to conclude the proof. However, for the convenience 
of the reader, we show here how to give a simple direct proof of the C 2,a regularity of u with 
a' = 2a. 

As in the previous step, it suffices to show that u is C 2,a at the origin, and for this we have to 
prove that there exists a sequence of paraboloids P k such that 

(5.23) sup \u-P k \ < CrJ (2+<T ' } 

B r§/C 

for some ro, C > 0. 

Let v k be as in the previous step, and let P k be their second order Taylor expansion at 0: 

P k (x) = v k {0) + Vv k (0) ■ x + ^D 2 v k {0)x ■ x. 
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We observe that, thanks to (|5.15|) . 

(5.24) \\v k - P h \\ CO{m ^— /R)) < \\ Vk - P k \\ c °(s hk+2 ) < C\\ D3 Vk\\c°(s hk+2 A /2 . 

In addition, by (I5.22D applied with j = 1, . . . , k and recalling that h k = h\2~ k and 2a < 1 (see 
(I5T20I) ). we get 

\\D 3 v k \\co [Shk+2 ) < ||D 3 vi||c°(S ha ) + H^ 3 ^ " D3v i+i\\c°(s hj+2 ) 

(5.25) k 

<c{i+j2 h{ r l/2) )< ch T i/2 - 



Combining pTl5]) . (f5T2l|) . (|5T25j) . and recalling (f57T7j) and ([5T20]) . we obtain 
so (|5.23p follows with tq = l/y/2 and a' = 2a. 



k ' 



□ 
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